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Abstract 

We devise a fast and optimal estimator for the amplitude of the bispectrum of clustered Infrared (IR) point-sources. We show how 
this estimator can account for the cases of partial sky coverage and inhomogeneous noise. Expected detection significance are pre- 
sented in terms of signal-to-noise, finding that the IR bispectrum will realistically be undetectable below 220 GHz with a Planck-like 
experiment; on the contrary detection may be achieved at, or above, 220 GHz especially if the CMB is removed. We also show how 
this estimator can be combined with estimators of radio and CMB non-Gaussianity to build up joint robust constraints. On the one 
hand, we find that, for a Planck-like experiment, CMB non-Gaussianity estimation can be decoupled from point-source contributions, 
unless few sources are masked. On the other hand, we find that the estimation of radio and IR non-Gaussianity are strongly coupled, 
which diminishes their separate detection significance. 
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1. Introduction 

The measurement and study of high order moments has emerged 
as an important field in observational cosmology as it probes 
deviations from Gaussianity. For example, observing non- 
Gaussianity (NG) from the Cosmic Microwave Background 
(CMB) or the large scale structure would probe primordial NG 
dating back to the generation of cosmological perturbations ; 
this in turn would constrain severely primordial models gener- 
ating these fluctuations (e.g. inflation). For example, a NG de- 
tection with the sensitivities of today's experiments would rule 
out a whole class of inflation models, among which the standard 
single-field slow roll inflation ( Acquaviva et al. 2003 ; Maldacena 



2003| lCreminelli & Zald arriaga 2004). 



The most popular form of primordial non-Gaussianity is the 
so-called 'local' type, parametrised by a factor /nl, for which 
the Bardeen potential takes the form : 



<D(x) = (D c (x) + / NL ((D G (x) 2 - <O c (x) 2 )) 



(1) 



where Op is the Gaussian part of the potential. It is on this 
form of NG that most efforts have focused and that the tight- 
est constraints were obtained. Indeed a fast estimators for / NL 
was developed and applied to the CMB data (Komatsu et al. 
|2005 2009, 201 1 1. However at CMB frequencies, many contam- 
inating signals are present, which typically have a non-Gaussian 
distribution. While the galactic emission is mainly confined to a 
defined area : the galactic plane, and may thus be masked, extra- 
galactic point-sources are present all over the sky and the faintest 
(below detection limits and confusion levels) are neither detected 
nor masked. These unresolved sources are therefore one poten- 
tially important source of secondary non-Gaussianity that needs 
to be estimated. 

There are principally two populations of point-sources of in- 
terest for the CMB : radio-loud sources with strong magnetic 
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fields yielding synchrotron or free-free emission in the radio do- 
main, and Dusty Star-Forming Galaxies (DSFGs) with thermal 
emission from dust heated by young star in the infrared (IR) do- 
main. Radio galaxies can be considered randomly distributed on 



the sky (Toffolatti etal.|1998 Gonzalez-Nuevo et al. 2005) at the 



CMB frequencies and thus have flat moments. Their skewness is 
related to their number counts as : 



bps 



dn 
dS 



■dS 



(2) 



where S is the source flux in Jy (bps may then be converted to 
CMB temperature elevation through Planck's law, but is usually 
quoted in dimensionless units AT/T in the literature). 
On the contrary, infrared galaxies are strongly clustered in dark 
matter halos, which typically enhances their correlation func- 
tions on large scales. Most recent CMB experiments, e.g. Planck 
([ Planck C ollaboration| |2006l l, ACT ( |Das et al.|[20TT] ) and SPT 
(Keisler et al. 201 1 ), are probing higher frequencies compared to 
previous missions so that the modeling and measurement of the 
infrared background becomes important. The Cosmic Infrared 
Background (CIB) and its anisotropics have indeed been de- 
tected and characterised at the power spectrum level at these 
frequencies (e.g. 



Hall et al.pOlO 
|20rib||Vieroetal.r2012D 



Planck Collaboration et al. 



Non-Gaussianity arising from point sources (radio or IR) is 
of particular importance for CMB studies, where maps are con- 
taminated by foregrounds which are imperfectly subtracted by 
component separation methods. Therefore estimation of the pri- 
mordial NG parameter /ml may be contaminated by the point- 
source NG ( |Babich & Pierpaoli| [2008} |Lacasa et aT[[20T2] i. 
Measurement of the non-Gaussianity of point-sources could also 
be a new way to further constrain the respective models for their 
emission and/or spatial distribution. 

A fast estimators for bps has been developed and applied to 
CMB data ( |Komatsu et al.|2005[[2009[[20TT) . Taking advantage 
of a prescription for the non-Gaussianity of IR source developed 
by ( |Lacasa et al.||2012"] l, the goal of this article is to propose a 
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fast estimator for the amplitude of the IR bispectrum. This can 
be used in non-Gaussianity studies to disentangle the primordial 
and foreground contributions so that measurements of both of 
them are unbiased. 

The article is organised as follows : in section [2] we briefly 
recall the prescription for the bispectrum of IR sources and in- 
troduce the estimator for its amplitude, in section [3] we tackle 
the problems of partial sky coverage and anisotropic noise, in 
section |4] we present the expected signal-to-noise ratio (SNR) 
of detection at several frequencies, in section [5] we show how 
to combine this estimator with the ones used for the CMB and 
radio galaxies to produce a joint constraint of non-Gaussianity, 
and we conclude in section [6] 

2. Bispectrum of IR sources 

2.1. Bispectrum 

Given a full-sky map of the temperature fluctuations A7\n) of 
some signal in direction n, it can be decomposed in the spherical 
harmonic basis 

a lm = j d 2 n r t Jn) AT(n) (3) 

with the usual orthonormal spherical harmonics Yg m 

J d 2 n T fffl (n) F* m ,(n) = 8 W <W (4) 

The power spectrum Q of the signal is the 2-point correlation 
function in harmonic space given as : 

(fltm a* e , m ,) - C( 5ec> 5 mm > (5) 

The bispectrum b^^j, is the 3-point correlation function in har- 
monic space: 

{a[ imi ae 2m2 a( }m ) - <J fl ^ 3 " x o^t^ (o) 
with the Gaunt coefficient 



n,2,3 



j d 2 n T 123 (n) 



(21 + l)i23 / e x h M/ d h h 
\n I II mi ni2 



(7) 



(8) 



where T, = Y(. m .. In the following the subscript 123 denotes 
the product of the corresponding variables e.g. X\r>, = X\ X2X3. 
G123 is zero unless the triplet (l\,ti, {3) follows the triangle in- 
equalities and mi + m% + = 0. 
The bispectrum estimator is : 



'123 



with 



— 7 G 



1,2,3 «1 «2 «3 



(2^+l)i23 / l\ h h 

~ 



(9) 



(10) 



being the 'number of modes' for the (t\,t% t l{) triangle. 

In the following, we note in bold the bispectrum which ac- 
counts for beam effect : 



We also note 



» 1,2,3 = b lxllh b (i bf 2 b ( 



C c = Ccb 2 c + C° oise 



The power spectrum accounting for both the noise and the beam 
effects. 



2.2. Optimal estimator of the amplitude of the bispectrum 

The non-Gaussianity from radio sources at CMB frequencies is 
entirely characterised by a single parameter bps which can be es- 
timated on a map ([Komatsu et al. 2005 ). Hence the radio-source 
contamination to /ml can be marginalised over. By contrast, the 
IR non-Gaussianity is more complicated and is characterised by 
the full bispectrum, whose estimation is computationally inten- 
sive (with O(N^) operations). 



In this context, Lacasa et al. ( 2012| l have proposed a full- 
sky analytical prescription to compute the bispectrum of extra- 
galactic sources from their number counts and power spectra. 
In this prescription, the bispectrum arising from a population of 
point sources writes : 



(11) 



with Cc being the power spectrum of the considered point 
sources, and 



3/2 



(12) 



where S is the source flux, the integral runs up to the detection 
limit of the resolved sources indicated by a flux cut S cut , and 4f 
are the number counts per steradian. 

This prescription gives a bispectrum in good agreement with the 
one measured on maps of radio and dusty star-forming galaxies 
simulated by Sehga l et aL] ( |2010| l. 

In the following, we will focus solely on IR point-sources. The 
level of non-Gaussianity induced by IR sources can thus be 
quantified by the amplitude of their bispectrum compared to the 
template bispectrum given by Eq[TT] The use of the amplitude 
of the bispectrum permits to derive easily the bias on primor- 
dial CMB estimation /nl, and to marginalise over this nuisance 
parameter as will be seen in Sect[5] 

The amplitude of the bispectrum is measured by minimising 
the x 2 of the map bispectrum to the theoretical one : 



with 



s 

t\<h<h 



(b 



1,2,3 



A h IR 



x A 



(13) 



(14) 



being the bispectrum variance in the weak non-Gaussianity 
limifM In this expression, Cf x is the map power spectrum, and 
we have 



16 equilateral triangle 
2 isosceles triangle 
1 general triangle 

Then the maximum likelihood estimator is: 

Zb 1,2,3 bj 2 3 „ 



(15) 



(16) 



1 This approximation is justified for the Cosmic Infrared Background 
(CIB) which is close to Gaussian (see 1 Lacasa et al. 2012) for details) 
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with 



3.2. Reducing variance and reaching optimality 



(T 



(]<(2<h 



(17) 



being a normalisation factor giving the error-bar in the assump- 
tion of a Gaussian likelihood, which is valid in the weak NG 
limit when <t 2 {(\, £2, £3) can be considered independent of A. 
The theoretical bispectrum, Eq[TT[ is separably ; hence a sim- 
pler computation of the numerator of Eq{l6]can be devised, in- 
spired by Komatsu et al.| ( |2005| ). Indeed, defining the filtered 
map: 



a 1 ' 3 Jcf-bt 
F(n) = ^ — a (m Y {m (n) 



/^>tOl 



and 
Sir 
we have: 



= j d 2 n F(n) 



b 1)2)3 6123 

tOt/^>tOt/^>tOt 



/ J /OtOt/^tOt/^tO' 

f,M 3 W 2 W 3 



x AT 123 = 6 x ^ 



bl l R 2,3 t>123 



(18) 



(19) 



(20) 



Furthermore it is worth noting that this approach is faster than a 
full bispectrum analysis, having O(N^) operations. 



3. Masked sky and inhomogeneous noise 

In realistic cases of CMB analysis, the statistical isotropy of 
the CMB signal is broken by inhomogeneous noise -e.g. due 
to scanning strategy- or by masking large areas of the sky -e.g. 
those contaminated by galactic emissions. In these cases, the es- 
timator of the IR-bispectrum amplitude is no longer optimal and 
it is biased in a non-trivial way. Nevertheless and similarly to the 
case of the /nl, the bias and the lack of optimality can both be 
tackled through adapted modifications to the estimator. 



3.1. Debiasing 

When a covariance matrix of the map(s) can be estimated, 
Creminelli et al. (2006) have shown that Wiener filtering of the 
map(s) will debias the non-Gaussian estimator from anisotropic 
contaminants or noise. Specifically, if C is the estimated covari- 
ance matrix in harmonic space, the estimator Eq. 16 can be de- 
biased by making the change 



a ln 



(C _1 -a), 

V Itm 



1 



at', 



(21) 



=-t..t 



in the filtered map Eq. [18] where a = (a( m )( = 2 
vector of harmonic coefficients. 

This step is sometimes called CT 1 prefiltering, and replaces 
the 1/Cf filtering of Eq"" 



18 



(in the isotropic case C 1 reduces to 
1 / Cf 1 ). An efficient al gorithm for this filtering has been designed 
by|Elsner & Wandelt|p012|l. 



2 a bispectrum is said separable if it can be written in the form 
b(( u ( 2 J3) = I li M£i)gi(t2)h i (£ 3 ) + perm. This is the case for the IR 
bispectrum with the sum reducing to 1 element and fj(£) = g t (€) = 



hi(i) 



,1/3 



r 1R 



When isotropy is broken, the 3-point correlation function used to 
define the bispectrum in Eq. |6]no longer has minimal variance. 
It must be replaced by the Wick product of the three harmonic 
coefficients given by: 

®i\m\ @?2tn2 ®£zm*± * ®t\m\ ^i^mi ®tzms - {ad 

) a (im .(22) 

This expression exhibits the same mean as in the isotropic case 
but it has a lower variance. Furthermore when applied to the 
isotropic case Eq. 22 gives the same bispectrum estimator as in 
Eq.|6](e.g. for l\ + ti + £3 all the expectation values vanish). 
Only the 2-point correlation function is to be considered for the 
expectation values. Therefore in practice, the latter are obtained 
from a sufficiently large number of Gaussian realisations with 
same power spectrum as that of the signal considered. 
Inputting the linear corrections given by Eq. 22 into the expres- 
sion of Sir in Eq. 20 yields: 



Sir -> S m = j d 2 n F(n) 3 -3 xj d 2 n F(n)G(n) 



with 



G(n) = JV /3 V C 'f C 'f < fl i fl 2>MC Ydn)Y 2 (n) 



= <F(n) 2 > MC 



(23) 

(24) 
(25) 



4. Expected detection significance 

In this section, we present the expected detection significance 
in terms of the signal-to-noise ratio for the previously described 
estimator. To this end, we use values of a and C' R that we mea- 
sure in the simulations of Sehgal et al. (2010) at frequencies 150, 
220, 280 and 350 GHz. The obtained values were shown to re- 
produce adequately the IR bispectrum at these frequencies (see 
Lacas aet al.| ( |2012"] l). 

4.1. Ideal cosmic variance-limited case 

Let us first focus on an ideal case where the signal is made solely 
of IR galaxies, i.e. noise and other astrophysical contributions 
are neglected. We will compute the signal-to-noise ratio of the 
detection of the non-Gaussian signal as defined by the amplitude 
of the bispectrum. Each triplet t\ < £2 < {3 contains information 
proportional to its number of configurations : 



^ 2 (A) 



l\<ti<h 



(26) 



In other terms, the SNR of the detection goes down with the total 
number of bispectrum configurations: 



is the SNR 



1 



cr(A) 



with 



N t . 



it(Anax) — / \ 



(27) 



(28) 



For example in the case of full-sky maps of the IR sources at fre- 
quencies 150, 220, 280 and 350 GHz with angular scales ranging 
from £ m i n = 2 to £ m . dx = 2048, the obtained SNRs are of order 
~1000 and they do not vary much from frequency to frequency 
(see first line of Table [TJ. 
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4.2. Realistic case 

In a real experiment, the IR signal/map is convolved by the in- 
strumental beam and the signal is contaminated by noise and 
other astrophysical signals, the most prominent being the CMB 
and the galaxy, which can be masked. Hence the SNR is signifi- 
cantly decreased compared to the ideal case. Table[TJsummarizes 
the expected SNR in the case of a 5 arcmin beam and 10~ 8 dT/T 
homogeneous detector noise across frequencies. We show how 
the detection level/significance decrease from the ideal case to 
the case of a convolved but "perfect" full-sky IR map, then in- 
cluding CMB contamination but still full-sky, and finally with 
a 50% sky-fraction mimicking the mask that may be applied to 
avoid galactic dust emission. We also include the case where 
90% of the CMB can be removed by a component separation 
method. For the partial-sky coverage cases, we assume that the 
optimisation and the debiasing described in Sect. 



been applied with the SNR scaling as /< 



-1/2 
SKY 



3.1 & 3.2 have 



frequency (GHz) 


150 


220 


280 


350 


Ideal case 


1218 


1157 


1161 


1159 


Full-sky IR 


15.5 


98 


336 


833 


Full-sky with CMB 


0.39 


6.7 


55 


387 


50% sky with CMB 


0.28 


4.7 


39 


274 


50% sky with 10% CMB 


3.67 


45.9 


210 


577 



Table 1. Expected SNR at £ max = 2048. Ideal stands for a full- 
sky cosmic-variance limited IR map without noise nor contami- 
nations. Full-sky IR represents the case where the IR map is con- 
volved with a 5 arcmin Gaussian beam and contains 10~ 8 AT/T 
noise (C[ = 7.4 • !0~ 4 fiK 2 ■ sr). The same instrumental specifi- 
cations are adopted for the following lines. Full-sky with CMB 
is when CMB is present, contaminating the IR map. 50% sky 
with CMB stands for the case of masking half of the sky. 50% 
sky with 10% CMB stands for a half sky IR map where 90% of 
the CMB (in amplitude) has been removed by some component 
separation. 



The SNR of the detection is severely reduced from ideal to 
more realistic cases especially at lower frequencies. The signal 
is still detectable at all frequencies when including beam and 
noise effects. However, when the CMB is included the signal 
becomes undetectable at 150 GHz. Above this frequency the de- 
tection significance is SNR~5<x at 220 GHz and reaches 40 to 
200cr at 280 and 350 GHz where the IR emission is dominant 
even in presence of CMB contribution. 

Removing the CMB, with a component separation method, per- 
mits to further improve the detection at 220 GHz by increasing 
the SNR to ~46<t. In this case, the IR non-Gaussianity at 150 
GHz seems marginally detectable. However at this frequency the 
non-Gaussianity from radio sources becomes an issue (while at 
higher frequencies IR sources are forecasted to dominate the un- 
resolved sources population). We tackle in the next section the 
issue of non-Gaussianity other than that of the IR sources. 



5. Joint NG estimation 

When several non-Gaussian signals are present, a joint estima- 
tion of their amplitudes taking their covariances into account 
is necessary. We focus in the present study on the main ex- 
tragalactic non-Gaussian signals, namely CMB, radio and IR 



point-sources. A joint estimation requires minimisation of the 
chi-square : 



X 



! (/nl>ps,A ir ) = 



(k 1*2,3 ~~ k™ d 3 el (/ N L,bpS,^IR)) 



(29) 



l\<t1<l% 
with 

b^ d 3 el (/NL, ^PS , Air) = / NL b™ B + b PS b^ + A IR bf 2j3 (30) 

and (involved in the definition of cr 2 {€\,{2, £3), see Eq|74"]> 

Cf = (cf MB + C RAD + Cf) b\ + C™ ise (31) 

We use the standard' parameter fops for compatibility with the 
literature (instead of noting e.g. Arad)- 

Then if we define the scalar product between two bispectra b" 
and b^ : 

^ y ^Agj (32) 

Minimising Eq|29]corresponds to solving the linear system : 



<b 



'CMB. "RAD 



> < b 



RAD' "RAD 



> <bi R ,b R 



\ 


f /NL ) 




bps 


/ 


< Air , 



( < bobs, b C MB > ' 
< b bs,bRAD > 

< b b s ,biR > ) 



(33) 



If we introduce -noted with upper tilde- the "naive" estimators 
which consider only one source of non-Gaussianity, e.g. 



iir 



< bpbsjbrR > 

< biR, biR > 



(34) 



then Eq(33]can be rewritten to define the Joint NG estimators 
-noted with upper hat- as : 



(35) 







'/>/ 


bps 


= AT 1 • 


fops 


, Air , 




k Air , 



with 
M = 
and 



1 



Afo™B 

4 



,IR 



AA§< B AA RAD 



1 

MR 



Ab 



PS 



(36) 



A fRAD _ < bRAD, bcMB > 
< "CMB? "CMB > 
< bcMB » b R AD > 



Ab 
AA 



CMB 



' S < bRAD, bRAD > 

CMB _ < bcMB, biR > 

< biR,biR > 



Ab 
AA 



NL ' 

IR _ 
PS 

! 

IR 



< biR, bcMB > 

< bcMB , bcMB > 

< biR, bRAp > 

< bRAD, bRAD > 

< bRAD, biR > 
< biR, biR > 



We compute for illustration the mixing matrix at 220 GHz. 
Neglecting the noise and beam effect and assuming the ERCSC 
flux cuts (Planck Collaboration et al. 2011a), we find for l m . dx - 
2048 : 



M = 



1 0.124 0.631 ) 
2.08-10~ 5 1 1.79 
3.16 10 5 0.535 1 



(37) 
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where we renormalised the radio amplitude to get numbers of 
order 1, i.e., the mixing matrix must be understood to be applied 
tofcps/(1.8-l(T 28 ). 

The matrix is mostly bloc diagonal, with the /ml part decoupling 
from the bps and Air parts. Therefore, unless bps » 1CT 28 , /nl 
estimation can be considered unaffected by extragalactic fore- 
grounds at this frequency. Conversely, estimation of NG from 
radio and IR is unaffected by /nl, but however we see that the 
two contributions affect heavily each other, with the £>ps-Atr sub- 
matrix being poorly conditioned (condition number ~ 180). 

When we include instrumental effects, namely a 5 arcmin 
Gaussian beam and a 10~ 8 Ar/T noise (same £ max and flux cut), 
we find the mixing matrix : 



M 



1 0.0461 0.245 ^ 
1.2110- 4 1 2.12 
1.31 10 4 0.435 1 



(38) 



The relative importance of CMB NG compared to point-sources 
NG has increased. Indeed beam and noise effects mostly de- 
crease the signal-to-noise in configurations with high multipoles 
; this down-weights point-sources NG which is dominant at 
high I, while CMB NG dominates at lower multipoles. However 
CMB and point-sources NG estimation can still be considered 
decoupled unless bps » 10~ 28 . In the £>ps-A ir submatrix, we 
see that the relative importance of bps has decreased compared 
to Air. Indeed the IR signal is most important at large angular 
scales while the radio SNR comes mostly from small angular 
scales which are affected by the beam smoothing and the noise. 

The fops-AiR submatrix is still ill-conditioned when we ac- 
count for beam and noise effect ; this is due to a high correlation 
between the IR and radio bispectra templates. This correlation 
can be quantified by the quantity r = cos 6 e [-1,1] : 



< b„,bfl > 



r afi 



(39) 



with a, ft = RAD,IR or CMB. In other terms Q afi is the angle be- 
tween the vectors b a , bp in the vector space of bispectra with the 
scalar product defined in Eq 32 At { max = 2048, with a 5 arcmin 
Gaussian beam and 10~ 8 Ar/r noise, the correlation coefficients 
of the different bispectra are listed in Table [2] 





IR/RAD 


IR/CMB 


RAD/CMB 


r 


96% 


0.57% 


0.24% 



Table 2. Correlation coefficients between the different bispec- 
trum templates at 220 GHz with £ max = 2048, 5 arcmin Gaussian 
beam and 10~ 8 Ar/r noise. 



In Fig{T] we plot these coefficients as a function of the maxi- 
mum multipole used in the analysis, with the same instrumental 
specifications. 

When analysing only large angular scales (small £ m - dx ), the 
radio and IR bispectrum templates are anti-correlated with the 
CMB bispectrum. This is because both radio and IR bispectra 
estimators are positive while the CMB is negative due to the 
Sachs-Wolfe effect term. However when going to higher resolu- 
tion, the first acoustic peak makes the CMB bispectrum positive 
with large values in some configurations (e.g. equilateral). This 
makes the correlation between the radio and CMB bispectrum 
positive for t > 168. The IR/CMB correlation in turn, becomes 
positive at slightly higher multipoles, i = 190, because the IR 
bispectrum peaks in the squeezed configurations where the CMB 




1500 2000 



Figure 1. Correlation coefficients between the radio, IR and 
CMB bispectra at 220 GHz with 5 arcmin Gaussian beam and 
10 -8 AT/T noise. Left panel : correlation between IR and CMB 
in red, between radio and CMB in blue, as a function of ( msx . 
Right panel : correlation between IR and radio. 



bispectrum is still negative. At larger multipoles, because of the 
complex pattern of the CMB bispectrum and the changes of sign 
with acoustic peaks, the correlation coefficients tend asymptoti- 
cally to zero. 

The IR and radio bispectrum templates are positively corre- 
lated as they are both positive ; this correlation is small when 
analysing only the lowest multipoles, but it increases rapidly 
over 80% with a small 'dip' at ^ max = 200 - 300 at the location 
of the first CMB acoustic peak. Indeed the CMB spectrum is in- 
volved in the denominator of Eq[32] so at £ m . dx = 200 - 300 it 
typically down-weights non-squeezed configurations where ra- 
dio and IR bispectra are the most correlated. Nevertheless, the 
correlation between radio and IR bispectra asymptotically tends 
to 1 at high angular resolutions, as the IR bispectrum flattens at 
high multipoles. 



8x10" 




Figure 2. Icr and 2cr confidence level for (Air , bps) in a joint 
estimation at 220 GHz £ mdx = 2048 with 5 arcmin Gaussian 
beam and 10~ 8 AT/T noise 



Finally, in Fig|2]we forecast the likelihood contours in the 
(Air , bps) plane, that we expect with the aforementioned instru- 
mental specifications. Air and bps estimations are quite degener- 
ate, as expected from the high correlations of their templates, 
which significantly degrades the constrain that can be put on 
each of them independently. However we can still have an or- 
der of magnitude estimation for both of them. Also, including 
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the prior bps > (and also to some extent A 1R > 0) would im- TofFolatti, L., Argueso Gomez, R, de Zotti, G., et al. 1998, MNRAS, 297, 117 

prove the constraints. Finally, CMB removal through component vier °. M - p -> w ™g. L -> Zemcov, M., et al. 2012, ArXiv e-prints 
separation may also help by decreasing the variance of both es- 
timators (although not their correlation). 



6. Conclusions and discussion 



The prescription proposed in Lacasa et al. ( 2012| l yields a sepa- 



rable form for the bispectrum of IR point-sources. This permits 
us to build a fast estimator for its amplitude a-la-KSW (Komatsu 



et al. 2005 1. We show how this estimator can be optimised to 
account for partial sky coverage and inhomogeneous noise. We 
compute the detection significance that can be expected at sev- 
eral frequencies. Finally we show how this estimator can be used 
together with estimators of the radio and CMB non-Gaussianity 
to build a joint estimation of sources of non-Gaussianity. We 
find that the estimation of CMB NG is barely correlated with 
the estimation of point-sources NG, so that unless the latter is 
too important -e.g. if few sources are masked-, both estimations 
can be considered decoupled. On the contrary, for point-sources, 
IR and radio NG have similar shapes especially at high mul- 
tipoles which make their amplitude estimation degenerate. We 
show how this hampers the detectability of the IR NG at 220 
GHz. Increasing this detectability can be achieved for example 
by subtraction of the CMB, or by including priors on bps in the 
analysis (e.g. positivity and/or upper limits given the flux cut and 
radio source models). 

This study opens up the possibility of robust measurements 
of non-Gaussianity where contaminations are treated in a sys- 
tematic way. It also opens up the possibility of the detection of 
point-sources NG, which may permit to constrain their respec- 
tive models. 
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